How robust are Majorana modes in multiband semiconductor wires? 
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We study the emergence of Majorana bound states (MBS) in multiband semiconductor wires with 
Rashba spin-orbit (SO) coupling, proximity coupled to an s-wave superconductor and in the presence 
of an external magnetic field. We consider long, but finite, multiband wires (namely, quasi- ID wires 
with dimensions Lx ^ Ly). Our numerically exact results demonstrate that interband mixing 
coming from Rashba SO (Rashba mixing) hybridizes Majorana pairs originating from different 
transversal modes while simultaneously closing the effective gap. Thus, regions were many MBS 
coexist are effectively trivial owing to Rashba mixing. On the contrary, Majorana physics is robust 
provided that only one single transversal mode, not necessarily the lowest one, contributes with a 
Majorana pair. We also find that Majorana physics is fragile with respect to orbital effects induced 
by any small out-of-plane component of the magnetic field. 

PACS numbers: 74.78. Na,71. 10. Pm,74.20.Rp, 74.75. -|-c 



I. INTRODUCTION 

Matter and anti-matter are related by a charge conju- 
gation operation whereby matter and anti-matter obey 
the relativistic Dirac equation with positive and negative 
energies, respectively^. Ettore Majorana proved that, 
under some symmetry considerations, the Dirac equa- 
tion could admit real solutions corresponding to particle- 
antiparticle pairs made by the same entity, a Majorana 
fermion^ . Majorana fermions were, firstly, proposed in 
the context of particle physics to describe neutrinos^ and 
more recently the Majorana search has been revived^^— in 
the condensed matter—"— and atomic physics^i^i^ com- 
munities. Apart from the fundamental interest of find- 
ing Majorana fermions, part of the excitement comes 
from the non-Abelian braiding statistics of these par- 
ticles which could be useful for topological quantum 
computation^ . 

The most important breakthrough in the condensed 
matter context originated with the seminal work of Fu 
and Kane who predicted that Majorana quasi-particles 
should appear in topological insulatorsi^ which, in prox- 
imity to s-wave superconductors, may behave as topo- 
logical superconductors. Similar ideas using semicon- 
ductors with strong spin-orbit (SO) couplingi^"— rapidly 
appeared. Arguably, the most promising semiconductor 
proposals are based on semiconductor nanowires (NW) 
with strong SO and giant Zeeman splittings, like InAs'^'^ 
and InSb^i^. This, together with the ability to induce 
superconductivity^^!^, has spurred a great deal of exper- 
imental activity towards detecting Majoranas in NWs. In 
one dimension (ID), these systems support a helical liq- 
uid state, where and external magnetic field perpendicu- 
lar to the SO field forces opposite spin states to move in 
opposite directions'^* . In the presence of superconduct- 
ing pairing, such helical state may become a topological 
superconductor supporting Majorana modes which ap- 




FIG. 1. (Color online) Sketch of a quasi-lD finite wire of 
dimensions {Lx, Ly) in close proximity with an s-wave super- 
conductor and in presence of a tilted magnetic field. 

pear as zero-energy Andreev bound states^^^. It can be 
shown^^ that this system is equivalent to Kitaev's model 
of ID spinless p-wave superconductors^. Such ideas have 
been recently generalized to the more realistic case of 
quasi-lD systems where the effect of transverse modes 
is taken into account. Both quasi-lD NWs^"— and 
p-wave superconductors^"— have been treated. From 
these works, an interesting picture emerges of successive 
filling of subbands giving rise to alternating trivial/non- 
trivial topological phases. The main idea behind this be- 
haviour is that multiband systems can support ID topo- 
logical superconducting phases with MBS provided that 
there is an odd number of partially occupied bands and 
that the width of the quasi-lD geometry does not ex- 
ceed the superconducting coherence length. Note, how- 
ever, that rigorous topology criteria can only be applied 
for infinitely long strips where Z2 invariants, such as the 
Majorana number—, unambigously determine the phase 
diagram. Thus, more finite-system studies beyond effec- 
tive two-band models^i^ and tight-binding^i^i^ are 
needed. Relevant questions such as the appearance of 
fermionic subgap states^ need to be addressed. Finite 
size calculations fully taking into account the SO inter- 
action are important inasmuch as the creation of long 
nanowires supporting a single propagating mode is still 
an experimental challenge. 
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We here address such task by performing numerically 
exact diagonalizations of the full SO model of a realistic 
quasi- ID nanowire in the clean limit. In particular, we 
study fundamental issues related to the robustness of Ma- 
jorana physics against {i) spin-orbit interband mixing, 
{a) finite size effects, and (Hi) magnetic orbital effects. 
Based on our results we draw two possible favorable sce- 
narios for the experimental detection of Majorana-like 
states. In the first one, Majorana physics is robust pro- 
vided that only one single transversal mode, not neces- 
sarily the lowest one, contributes with a single Majorana 
pair. This condition originates from interband Rashba 
mixing which strongly hybridizes the Majorana modes 
in situations where more than one pair of them coex- 
ist. Alternatively, we show that systems with negligi- 
ble SO mixing may host many zero energy modes. This 
second scenario can be achieved in NWs without inter- 
band SO interaction. Such situation can be engineered in 
NWs without SO where either curvature or an inhomo- 
geneous Zeeman field generate an effective band-diagonal 
SO couplin^^. We also stress the importance of the ori- 
entation of the magnetic field for the Majorana physics 
detection. Any misalignment with the plane of the NW 
introduces undesirable magnetic orbital effects that, as 
shown below, affect dramatically the formation of Majo- 
rana modes. 



II. QUANTUM WIRE MODELIZATION 

We consider a quasi- ID semiconductor strip with a 
strong Rashba SO and proximity-induced pairing, in 
the presence of a magnetic field B arbitrarily oriented 
(see Fig. [J). The model Hamiltonian can be written as 
Hnw = y-sp + T-Lz + T-Lso-, where 

IP +IP 

nsp^ ' ^^i + V{x;La,) + V{y;Ly) , (1) 
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"Hz = -gsMfisB ■ a , 
Uso = ^ (? X n) • f 



(2) 
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with TO* the effective mass, Hi — pi + {e/c)Ai the canon- 
ical momentum (e is the electrons charge, and c is the 
light speed) and A the potential vector. Hsp is the 
single particle energy, refered to the chemical poten- 
tial /i, in a hard-wall confinement in both x and y di- 
rections, given by V{a; L) = for < a < L and 
V{a] L) = oo otherwise. Hz is the standard Zeeman 
term, with = qsmI^bB 12 being the Zeeman energy. 
T-Lso is the Rashba SO Hamiltonian. The constant vec- 
tor 8 in T-Lso corresponds to the effective electric field 
due to confinement. In our quasi- ID geometry, typical of 
Rashba systems, the strongest confinement occurs along 
the z axis and, therefore, Hso ~ oi (a^IIj, — Oylix) with 
a = In a convenient gauge, A = —Bzyux, it is 

IIx = Px — &yBz/c and Ily = py. Notice then that both 



"Hsp a-nd T-Lso a-re modified by the magnetic orbital ef- 
fects, represented by the magnetic length — \/fic/eBz- 
More specifically, in the case of the Rashba Hamiltonian 
we have 
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so 
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where the distinction between ax and ay is artificially 
introduced. Clearly, in Rashba NWs how- 
ever in a recent proposal situations with ay ^ ax, includ- 
ing the case ay = 0, can also be engineered^. This is 
achieved in NWs without SO interaction using inhomoge- 
neous fields B{x) or modulated g-factors gsAiix). In gen- 
eral, Eq. Q mixes different quantum well subbands in x 
(through the px operator) and in y {py and y operators). 
Assuming the condition Lx 3> Ly the relevant mixing 
term is UyPyGx that we name Rashba mixing. Finally, we 
stress that orbital effects modify the SO coupling, includ- 
ing the Rashba mixing, and, as shown below, alter dra- 
matically the topological phases of multiband nanowires. 
The second-quantized Hamiltonian can be written as 



Ti 



NW 



where n = 

{nx,ny} are the quantum numbers associated with the 
transversal modes due to confinement and a =t, \. de- 
notes the spin. When the nanowire is proximity-coupled 
to an ordinary s-wave superconductor we add the BCS- 

and obtain 



Hamiltonian, Tisc — Ylr 
the low-energy physics of this system in a Bogoliubov 
deGennes description T-LbcIG ~ 5 
where '^n = (ci-^., c,\;, c„;, -c„t-)1' and 
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Our results below, obtained by exact numerical diago- 
nalization of Eq. ([S|), demonstrate the emergence of NW 
gapped low-energy eigenstates; ideally zero modes sepa- 
rated by a sizeable energy gap from the rest of the spec- 
trum and, as anticipated, we will focuss on the robustness 
of this behavior. 



III. MAJORANA-LIKE BEHAVIOR IN FINITE 
NANOWIRES 

As good estimates, let us notice that in infinite ID 
nanowires the occurrence of a Majorana pair requires 
large enough Zeeman energies with fields exceeding a 
critical value Be given hy gsMfJ-BBc/2 = \J [i^ + A ^-^^'-^^ . 
In the quasi-lD case, a similar condition for the critical 

field Be " of the Majorana pair from each Uy transverse 

mode is gsuP^BB'^J'''^ 12 = V(£«„ - [if + A^, with e„^ 
the transverse mode energy. Our results of Fig. 2, corre- 
sponding to the absence of Rashba mixing, confirm this 
scenario in long-enough wires and for field along x. We 
clearly see that increasing the field strength, successive 
gapped low-energy pairs, each belonging to a different 
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FIG. 2. (Color online) Results neglecting Rashba mixing, 
Oy — in Eq. (Q, and without magnetic orbital effects, 
for field orientation along x. (a) Spectrum of eigenvalues 
as a function of Zeeman energy As = gsMf^sB /2. Only 
the twelve eigenvalues closer to zero are displayed; higher or 
lower eigenvalues than the ones displayed correspond to bulk 
excitations and are not shown. Parameters: Ly — 150 nm, 
L:,, = 3 ^m, A = 0.075 meV, n = 0, = 0.045 eVnm. 
(b) Probability density for the Majorana-like edge modes cor- 
responding to the first and second transversal modes for a 
Zeeman energy of 2 meV in panel a) . 



transverse mode, appear in the spectrum. This picture 
is consistent with the results of Refs42ri2. We empha- 
size, however, that Fig. 2 corresponds to a situation in 
which transverse modes are uncoupled. Moreover, the 
pair energies are not exactly zero, but oscillate around 
zero with increasing B, showing a steadily increasing 
amplitude. This is a finite size effect, indicating that 
the shortness of the wire eventually dominates, remov- 
ing the pairs from zero energy and destroying the gap 
with the nearby states as well. The existence of gapped 
low-energy modes (ideally zero modes) and particle-hole 
symmetry are the requirements for the emergence of Ma- 
jorana physics. Let us label positive energy states at a 
given B in Fig. 2 by an index / = 1,2,... in increas- 
ing energy order. Analogously, negative energy states in 
decreasing energy order are labelled by / = —1,-2,.... 
According to Dirac, positive-energy states are particle 
states while negative ones are their conjugated or an- 
tiparticle ones. Clearly, when \Ei\ is sizeably nonzero 
|/) and I — /) are different stationary eigenstates of H, 
particle-antiparticle conjugates of each other. Notice now 
that we may form the two combinations 



7(^)) = ^(|/) + |-/)), 
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with the remarkable property that each of them is its own 
antiparticle conjugate. These are the Majorana states 
that, in general, are not Hamiltonian eigenstates. How- 
ever, a conspicuous exception occurs when Ej « E^i « 

since then 17^^'') and |7^^'') are approximate zero eigen- 
states of v., as occurs with E±i when > 0.5 meV 
in Fig. [2^. In this case the pair of Majoranas are gap- 
protected quasieigenstates and thus robust and almost 
stationary in time, what makes them interesting from 
the applications point of view^. 



IV. ROBUSTNESS AGAINST RASHBA MIXING 

As we discussed in the previous section, in the ab- 
sence of Rashba mixing, ay in Eq. ^ , successive gapped 
low-energy pairs, each belonging to a different trans- 
verse mode, appear in the spectrum. As we dis- 
cuss now, Rashba mixing makes this so-called even-odd 
behaviour—^— fragile when many transverse modes con- 
tribute to Majorana formation in realistic finite samples 
{LjLy = 20 and Ly « 4^50 ~ 0.4^, with Iso = Sli 
and e = ~ m'ltiso being the SO length and the 

superconducting coherence length, respectively). Panels 
a and b in Fig. [3] summarize the evolution of a repre- 
sentative spectrum when the mixing is gradually intro- 
duced. This mechanism leads to an effective coupling be- 
tween low energy modes, with relevant modifications of 
the spectra in the regions that for Fig. [2] contained two 
and three Majoranas. After the collapse of the second 
mode, a small gap near zero energy appears in Fig. [3^,b. 
More importantly, the presence of Rashba mixing tends 
to close the separation gap betwen low energy and higher 
energy states. These effects are thus quenching the ro- 
bustness of the Majorana modes by making them loose 
both their zero-enery character and their gap-protection 
from the rest of the states. Overall, all the regions with 
more than one Majorana pair, namely > 0.75 meV 
for Fig. [3Jd, become trivial. The density distributions of 
the three lowest eigenstates for As —4 meV is shown 
in Fig. proving that the character of edge modes is 
lost in presence of full Rashba mixing. On the other 
hand, the region with one Majorana pair, approximately 
in the range 0.25 meV < A^ < 0.75 meV in Fig. 
is clearly gapped. We have performed detailed calcula- 
tions as a function of the chemical potential (Appendix) 
that prove that these regions with a single Majorana 
pair are robust. Moreover, they show alternating triv- 
ial and non-trivial behavior as the chemical potential 
increases^2ii^d&. This alternating behavior appears as 
successive transversal modes are filled and a single Ma- 
jorana pair, originated from the highest mode, appears. 
Therefore: 1) Majoranas do not necessarily emerge from 
the first transverse mode and 2) they are robust even 
in regions Ab, namely in situations where a single 
band NW would be trivial. Our results place the Rashba 
mixing as ultimate responsible of the absence of clear 
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FIG. 4. (Color online) a) Energy of the lowest eigenvalue as a 
function of the field tilting angle (circles) and of the gap with 
the next state (squares). We have used a Zeeman energy of 
0.45 meV. b) Same as Fig.[2]for a field tilting angle of 1 degree 
out of the plane. Same parameters of Fig. [3)d. 
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have considered, we do not find such alternating behav- 
ior between trivial and non-trivial phases by varying the 
Zeeman field, as opposed to Ref.— . Nevertheless, even- 
odd physics in the whole (^, Ab) phase spaced"— should 
be recovered for long enough samples where topological 
arguments become exact^". We have checked that this is 
indeed the case by performing Pffafian calculations (Ap- 
pendix). Note that, of course, the contrary is not true 
and PfafSan calculations, valid when momentum along 
the wire is a good quantum number, do show non-trivial 
regions for values of fi and where the exact calcula- 
tions for realistic finite-size wires predict trivial behavior. 



FIG. 3. (Color online) Same as Fig. [2] including Rashba 
mixing. The same parameters have been used, except for 
ay = O.Gaa; in panel (a) and ay — a^ in panel (b). The 
densities of panel (c) are for the three lowest modes when 
As = 4 meV. 



sequences of nontrivial/trivial phases for large magnetic 
fields. Thus, a good strategy towards detecting Majo- 
ranas is to measure at magnetic fields close to the criti- 
cal field when the first transversal mode contributes with 
one Majorana pair. In such situation, the successive fill- 
ing of high transversal modes, as one increases the chem- 
ical potential, will produce alternating nontrivial regions 
with a single Majorana pair each, so that Rashba mixing 
is not efficent. We emphasize that, for the lengths we 



V. ROBUSTNESS AGAINST ORBITAL 
EFFECTS 

All the above results were obtained assuming a mag- 
netic field with a perfect in-plane orientation. We finish 
this discussion by studying the role of magnetic orbital ef- 
fects, due to Bz 7^ 0. The polar angle 6c = tan~^ Bz/ 
quantifies an unwanted small out-of-plane component of 
the magnetic field. Our results are shown in Fig. 2] A 
rather small vertical component of the field is sufficient 
to steadily increase the energy of the lowest mode, from 
a clear Majorana-like character with i?i « for Oc = 0° 
to El K, 0.2 A for Oc — 3°. At the same time, the gap 
from the lowest to the next state, E2 — Ei, is reduced 
as the magnetic field inclination decreases. Thus, orbital 
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effects severely suppress the Majorana character of the 
low-energy modes. This is demonstrated in Fig. Hl^b), 
where just one degree of tilting is enough for the spec- 
trum not to show gapped low energy modes as a function 
of B. 



VI. CONCLUSIONS 

We report the existence of Majorana modes in quasi- 
ID semiconductor wires. We draw two possible practi- 
cal scenarios for the detection of this gapped low-energy 
eigenstates: (i) Rashba NWs hosting a single Majorana 
pair; (ii) NWs with negligible SO interband mixing which 
can host many zero energy modes. This second sce- 
nario can be engineered in NWs without interband SO 
interaction^. Finally, we have addressed the fragility of 
the Majorana behavior against magnetic orbital effects, 
finding that a small misalignment of the magnetic field 
with respect to the NW is sufficient to quench Majorana 
physics. 
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FIG. 5. (Color online) Spectrum of eigenvalues as a function 
of chemical potential /i. Upper panel corresponds to the case 
without Rashba mixing , Uy — 0: (a) As = <?SAf/isB/2 = 
0.48 meV, and (b) As = 2.4meV. Lower panel shows the 
case with Rashba mixing, = a^,: (c) As = 0.48meV, 
and (d) As = 2.4meV. Only the twelve eigenvalues closer to 
zero are displayed. Parameters: Ly = 150 nm, = 3 /im, 
A = 0.075 meV, = 45 meVnm. 
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Appendix A: Supplementary Information 

1. Robustness of Majorana-like states against 
chemical potential 

In this appendix we present results for the nanowire 
spectrum as a function of the chemical potential. These 
results demonstrate that, in the presence of Rahsba mix- 
ing {ayPyUx term in Eq. 4 of the main text), regions with 
a single Majorana pair are robust while regions where 
Majorana pairs coexist are not. To illustrate this, we 
first show the spectrum without Rashba mixing (Fig. 5, 
upper panel). This energy spectrum can be understood 
as a superposition of the partial spectra corresponding 
to subspaces of the transversal modes Uy = 1,2,3- ••. 
For low B fields (Fig. 5a for As = 0.48meV) single 
Majorana pairs never coexist: as /i increases, successive 
traversal modes ny = 1,2,3- •• become populated and 
contribute with a Majorana mode each. These regions 
with Majoranas are separated by gapped regions without 
zero modes. These intermediate regions correspond to so- 
lutions with contributions from different Uy (for instance, 
the region around ^ — 2.5 meV corresponds to solutions 
where both Uy = 2 and Uy = 3 have a significant weight). 
At high magnetic fields (Fig. 5b for = 2.4meV) 
Majorana-like states from different transversal modes 



(i.e., different Uy) have a substantial overlap. One, two 
and three Majoranas-like states can coexist as shown in 
Fig. 5b (the region with three Majoranas is distorted 
because of finite size effects). We emphasize, again, that 
this coexistence of Majorana modes is possible owing to 
the absence of SO mixing. As we discuss below, a full 
calculation (Fig. 5d) shows that these regions are indeed 
trivial. For low B fields, where only one Majorana pair 
appears in the spectrum, the physics is similar without 
or with Rashba mixing (Fig. 5c for As = 0.48meV). 
Namely, alternating trivial/nontrivial regions separated 
by gapped regions without zero modes. Interestingly, 
non-trivial regions occur even for large chemical poten- 
tials /i ^ Ab (e. g. the region around /i = 3 meV in 
Fig. 5c). This result is in agreement with Refs^S"— 
. However, at high B fields (Fig. 5d for Ab = 2.4 
meV), the coexistence of various Majorana states is re- 
moved because of the Rashba mixing. All modes strongly 
hybridize, gapped regions disappear and the full spec- 
trum is trivial. We thus conclude that in the presence 
of Rashba coupling, the only robust Majorana regions are 
those with single Majorana pairs. Therefore, alternat- 
ing trivial and non- trivial regions can be expected in 
situations where only a single Majorana pair exists. Im- 
portantly: 1 ) this Majorana does not necessarily emerge 
from the first transverse mode and 2) they are robust 
even in regions /i 3> Ag, namely in situations where a 
single band nanowire would be trivial. These results sup- 
port the fact that Rashba mixing is the ultimate respon- 
sible of the absence of clear sequences of nontrivial/trivial 
phases for large magnetic fields. 
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FIG. 6. (Color online) Phase diagram (Pfaffian) for an infinite 
calculation versus As — g^BB/2 and chemical potential /i. 
Dark regions indicate nontrivial topological phases whereas 
light ones correspond to trivial phases. Parameters: Ly — 
150 nm, A = 0.075 meV, = ay = 45 meVnm. 



2. Phase diagram in the Lx ^ oo limit 



Our results are obtained by finite-size numerically ex- 
act calculations so it is interesting to investigate to what 
extent they agree with exact topological phase diagrams 
in the ^ oo limit. In order to do so, we have also 
performed Pfaffian calculations^ for an infinitely long 
nanowire (Fig. 6). Using the Pfaffian one predicts al- 
ternating trivial/non-trivial regions in the phase space 
(/ijAs). At low fields, this even-odd behavioiiiS"— agrees 
with our full results where the filling of successive Uy 
contributing with a Majorana mode each gives robust 
non-trivial regions, even with SO mixing. At high fields, 
however, the Pfaffian still shows even-odd behavior in 
regions where the spectra of our full calculations gives 
strongly hybridized modes and, hence, trivial behaviour. 
Thus, Pffafian calculations seem to become unreliable in 
regions where many transversal modes contribute with 
Majoranas and overlap. This can be understood as finite- 
size effects, very relevant in the case of strong SO mixing, 
are completely neglected. In such situations, full calcu- 
lations, as the one presented here, are needed. 
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